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We present a proposed design for a pair-breaking photodetector for far-infrared and sub-minimeter 
radiation. Antenna-coupled radiation generates quasiparticles in a superconducting absorber, the 
density of which are measured using a single Cooper-pair box. Readout is performed using an 
electromagnetic oscillator or a microwave resonator, which is well suited for frequency multiplexing 
in large arrays. Theoretical limits to detector sensitivity are discussed and modeled, with predicted 
sensitivities on the order of 10~^^ W/a/Hz. We anticipate that this detector can be used to address 
key scientific goals in far-infrared and sub-millimeter astronomy. 



Cryogenic radiation detectors based on low- 
temperature superconductors have a long and successful 
history in a wide variety of applications, including 
astronomy,^ high-energy physics,^ optical science,^ and 
communications^ over a spectral range from x-rays to 
the sub-millimeter. Various physical processes have 
been exploited to create sensitive detectors, which can 
be roughly separated into bolometric devices, such 
as the transition-edge sensoilSl and the hot-electron 
bolometer mixer and pair-breaking devices, such as 
the superconducting tunnel junction (STJ) detecto]!^ 
and the microwave kinetic inductance detector (MKID) 
P In a pair-breaking detector, radiation coupled to a 
superconducting absorber breaks Cooper pairs in the 
material, generating quasiparticles. A key challenge in 
designing pair-breaking detectors is the readout, the 
means of measuring the number of quasiparticles in 
the absorber in the presence of the Cooper pairs. This 
has been accomplished by measuring a current across 
a tunnel junction, as in the STJ, and by measuring a 
change in the impedance of a superconducting resonator, 
as in the MKID. 



In this work, we propose a pair-breaking detector ap- 
propriate for sub-millimeter radiation based on the sin- 
gle Cooper-pair box (SCB), a mesoscopic superconduct- 
ing circuit which has been extensively studied in the 
context of quantum computation.^ The SCB has also 
been shown to be extremely se nsitive to the presence 
of non-equilibrium quasiparticles.HE^E^^The proposed 
device, which we call the quantum capacitance detector 
(QCD), uses the SCB to sample the density of quasipar- 
ticles in the absorber. The QCD promises excellent sen- 
sitivity, with estimates of the minimum noise-equivalent 
power (NEP) on the order of 10~^^ W/a/Hz under real- 
istic operating conditions. Like the MKID, the QCD can 
be read out in the frequency domain, naturally lending 
itself to large-scale multiplexing using well-established 
fabrication and measurement techniques. Among many 
other potential applications, large arrays of sensitive low- 
noise detectors are critically needed to meet scientific 
goals in far-infrared and sub-millimeter astronomy.!^ 
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FIG. 1: (Color online) A) Cartoon schematic of the QCD 
concept. Incident submillimeter radiation is coupled by an 
antenna into a superconducting absorber, generating quasi- 
particles. The quasiparticle density in the absorber is then 
measured with an SCB. B) Circuit diagram of the SCB and 
its readout. 



I. DETECTION CONCEPT 

A. Overview 

The SCB consists of a small superconducting island 
coupled to superconducting leads (hereafter referred to 
as the absorber) by a pair of ultra-small Josephson junc- 
tions, typically arranged in a DC-SQUID configuration. 
A schematic of the SCB is shown in figure 1. The ca- 
pacitance of the island, which is typically on the order 
of a few fF, is small enough that the charging energy 
is the dominant term in the Hamiltonian, and charge 
number states are well defined. As the potential of the 
island is adjusted with a gate capacitor, the two low- 
est energy levels form an avoided level crossing. At the 
charge degeneracy point, where the lowest-lying eigen- 
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states are symmetric and antisymmetric superpositions 
of charge states, the SCB may be treated as a two-level 
quantum system coupled to an effective magnetic field. 
In the QCD concept, the state of the SCB is read out 
dispersively using a quantum capacitance technique, l^^^^ 
where the SCB island is coupled capacitively to an LC os- 
cillator which is slow compared to the qubit energy level 
splitting. By probing the oscillator with RF refiectome- 
try, one can directly measure the quantum capacitance of 
the SCB, which is proportional to the second derivative of 
the SCB energy level with respect to the gate charge. An 
equivalent measurement can also be performed using an 
inductive coupling between the SCB and the oscillator.^ 
Note that in the QCD scheme, the SCB is always oper- 
ated in its ground state, obviating the need for coherent 
manipulation. 

Although experiments aimed toward quantum compu- 
tation involve delicate superpositions of Cooper pairs, the 
SCB is exquisitely sensitive to incoherent quasiparticle 
tunneling between the absorber and the island, even at 
temperatures where thermally excited quasiparticles are 
strongly frozen out. When non-equilibrium quasiparti- 
cles tunnel across the SCB junctions, the parity of the 
device abruptly switches, destroying quantum coherence. 
As such, this effect is typically known as quasiparticle 
"poisoning" . However, this property can be exploited to 
make an extremely sensitive measurement of the quasi- 
particle density in the absorber, since the quasiparticle 
tunneling rates and the fraction of time the device spends 
in each parity state are related in an extremely simple 
way to the quasiparticle density in the absorber. 

The architecture of the QCD is shown in figure 1. 
Infrared or sub-millimeter radiation coupled to the an- 
tenna provides energy to break Cooper pairs in the ab- 
sorber, generating phonons and nonequilibrium quasipar- 
ticles which diffuse to the junctions of the SCB. In alu- 
minum, an attractive material for fabrication of SCBs 
and superconducting resonators, pair-breaking can occur 
for photons with a frequency above 90 GHz. The SCB 
is capacitively coupled to an LC oscillator or microwave 
resonator, which is used to measure the SCB state. 

When a single quasiparticle tunnels onto the SCB is- 
land, it switches the parity, changing the effective gate 
charge by one electron, bringing the SCB far from its de- 
generacy point and changing its capacitance. This sud- 
den shift in the capacitance leads to a large frequency 
shift in the oscillator, which was measured to be 140° 
in a recent experiment,'^ and is shown in Fig. 2. When 
quasiparticle tunnels back into the absorber, the parity 
switches back and the device returns to the degeneracy 
point. This results in a time trace given by a two-rate 
random telegraph signal. By measuring the average cen- 
ter frequency of the oscillator, one can average over this 
telegraph signal and extract the fraction of the time the 
SCB spends in the "even" and "odd" parity states, which 
is simply related to the tunneling rates and hence the 
density of quasiparticles in the absorber. In this paper, 
we refer to the even parity state as the one with zero 



^ 100| ■ ' • 



B) 



1000 



.x10 



2000 3000 
Time (|as) 



4000 



5000 




-150 



c) 



-100 -50 50 100 
Oscillator Phase (degrees) 

Log Number of Counts 



150 




Gate Charge (electrons) 



FIG. 2: (Color online) A ) Example of two-rate telegraph noise 
in the capacitance signal, recorded at the SCB degeneracy 
point. Data adapted from Ref. (14). B) Phase shift his- 
togram, taken at the degeneracy point, showing two clearly 
separated peaks corresponding to odd and even parity states. 
Note the large (140°) phase shift between the two. The phase 
shift signal in the QCD is an average taken over the entire time 
signal, as indicated by the red line. C) Phase shift histograms 
as a function of gate voltage. The z-axis is plotted logarithmi- 
cally to emphasize the presence of the minority (odd) parity 
state. The white lines are theoretically predicted capacitance 
traces for the odd and even parity states. 



or an even number of quasiparticles on the island, and 
the odd state as the one with one or an odd number of 
quasiparticles on the island. 

In the QCD concept, the absorber is small compared 
to the characteristic diffusion length of quasiparticles in 
Al, so that the spatial distribution of quasiparticles in 
the absorber may be taken to be uniform, and any given 
quasiparticle may be assumed to rapidly sample the en- 
tire volume of the absorber. Quasiparticles generated in 
the higher-Tc antenna will be trapped in the absorber 
by Andreev reflection. A key advantage of the SCB as a 
detector readout is that it is sensitive to the density of 
quasiparticles, rather than the overall number. Thus the 
absorber volume is a key design parameter which can be 
used to tune the operating wavelength range, sensitivity, 
and saturation power, as discussed below. 

The QCD concept is related to but clearly distinct from 
previous proposals to use single-charge devices for elec- 
tromagnetic detection. In one such experiment involv- 
ing a single-electron transistor (SET) with a supercon- 
ducting island and normal leads, Andreev-cycle transport 
was switched on and off when a quasiparticle from the 
leads was transferred to the island via photon-assisted 
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tunneling .EHl In a second experiment, quasiparticle tun- 
neling current in a biased submicron SIS junction is mea- 
sured using an RF-SET.^^ However, detectors based on 
the superconducting SET present significant challenges 
for multiplexed readout, due to the need to transform 
the high impedance of the SET junctions to that of a 
50 Q transmission line. 



B. Strategies for Multiplexed Arrays 



Another significant advantage of the QCD concept is 
the natural use of frequency multiplexing to read many 
pixels simultaneously with a single high-frequency line. 
Each pixel in a QCD array would have a different reso- 
nant frequency, determined by the values of the on-chip 
inductor and capacitor or the length of the resonator. 
All pixels could be read out simultaneously by applying 
a frequency comb to all devices through a common trans- 
mission line. Each frequency in the comb would match 
that of a single resonant circuit, and the refiected power 
at that frequency would constitute the SCB readout sig- 
nal for that pixel. Signal generation, demodulation, and 
analysis for large arrays can be performed using technol- 
ogy developed for software defined radio, and such mea- 
surement techniques are already in development for mul- 
tiplexed MKID arrays. A design for the electronics re- 
quired to implement such a multiplexing scheme is shown 
in Fig. 3. A simple multiplexed quantum capacitance 
readout of two SCBs has been recently demonstrated.^^ 

Since the operating point of the SCB is sensitive to 
small static charge inhomogeneities, it is impractical to 
tune the operating points of all SCB pixels in a QCD 
array to the degeneracy point with a single DC gate line. 
As shown in Fig. 2, the response of the SCB to parity fiuc- 
tuations is maximal at the degeneracy point. However, 
this problem can be readily solved by averaging over all 
gate voltages, as shown schematically in Fig. 3. An AC 
voltage with an amplitude of e/Cg (where Cg is the gate 
capacitance) can be applied to all SCB gates through 
a common line. This adiabatically sweeps the operat- 
ing point of each SCB through one full period, modu- 
lating the RF output signal. A mixer demodulates the 
refiected RF signal at the AC modulation frequency, and 
a down- converter translates the result to DC. In this way, 
a single AC tone applied to all SCB gates will appear as 
a sideband of each RF readout frequency in the comb, 
which gives the time-averaged capacitance averaged over 
all gate voltages when mixed down to DC. The refiected 
RF comb, containing the phase shift information for the 
entire array, can be demodulated at the gate modulation 
frequency, down-converted to the 0-200 MHz band, then 
digitized and digitally demultiplexed. 
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FIG. 3: (Color online) Schematic block diagram of proposed 
QCD readout electronics. Inset, lower right: Single pixel - An 
RF carrier signal excites the tank circuit of a particular detec- 
tor pixel, producing a reflected signal that is phase shifted ac- 
cording the average charge state of the SCB. An AC bias volt- 
age is applied, with an amplitude of e/Cg, which modulates 
the phase shift signal by sweeping the SCB through one full 
gate period. A mixer demodulates the reflected RF signal at 
the bias voltage modulation frequency, and a down- converter 
translates the result to DC. Main: Multiplexed readout - A 
low frequency comb function (0-200 MHz) containing several 
frequency components is produced digitally and block up- 
converted, resulting in a comb of RF carrier frequencies, each 
corresponding to a particular detector. All of the SCB gate 
lines are tied together and modulated at the same frequency. 



II. DEVICE PHYSICS 

A. Quasiparticle Tunneling 

The physics of nonequilibrium quasiparticle tunneling 
in single Cooper-pair devices has been extensively stud- 
ied experimentally and can be effectively described by a 
kinetic trapping theory.'^ With the SCB biased at the 
degeneracy point, the island will behave as a potential 
well or "trap" for nonequilibrium quasiparticles in the 
absorber, with an effective depth 6E = Ec — ^ — A, 

2 

where Eq — 2^ is the single-electron charging energy, 
Cy. is the overall capacitance of the SCB island, Ej is 
the Josephson energy, which is proportional to the crit- 
ical current of the junctions, and A = A/ — A^ is the 
island-absorber superconducting gap profile. While in 
superconducting circuits for quantum computation A is 
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frequently engineered to be as large as possible to sup- 
press quasiparticle tunnelingj^ in the QCD it is advan- 
tageous to have A < 0. For the remainder of this paper, 
we assume that A = 0. 

Given a uniform nonequilibrium quasiparticle density 
riqp in the leads of the SCB, quasiparticles will tunnel 
onto the island at a proportional rate Fin = Kuqp. The 
proportionality constant is given by 
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where Nl = D{EF)V^7TAkBT is the density of quasipar- 
ticle states available in the lead, D{Ef) is the normal- 
metal density of states at the Fermi level, Gat is the tun- 
neling conductance, T is the sample temperature, and 
/cb is the Boltzmann constant. Furthermore, the effec- 
tive density of states is given by 
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which includes BCS coherence factors accounting for 
quantum interference between electron-like and hole-like 
quasiparticle tunneling. The approximation of a tun- 
neling rate proportional to the quasiparticle density is 
valid as long as the operation temperature /c^T <C A. 

To a first approximation, the tunnel rate Fout for 
nonequilibrium quasiparticles on the SCB island back 
onto the absorber is independent of the quasiparticle den- 
sity riqp. The probability that the SCB will be found in 
the odd parity state then depends on the quasiparticle 
density in a particularly simple way. 
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A theoretical expression for Fout , which involves two sep- 
arate time scales for elastic and inelastic quasiparticle 
tunneling, can be found in Ref. (|^. In the proposed 
detector, Podd is the experimentally accessible quantity 
which is used to measure the quasiparticle density. At 
100 mK, the proposed QCD operating temperature, both 
K and Fout are only weakly dependent on temperature.^^ 
They are primarily set by the material and qubit param- 
eters. 



B. Oscillator Response 

When a quasiparticle generated by incident radiation 
tunnels onto the SCB island, the SCB parity switches, 
changing the quantum capacitance Cq. This shift in 
capacitance leads to a shift in the center frequency of 
the oscillator to which the SCB is capacitively coupled. 
Within the bandwidth of the oscillator, these shifts can 



be observed directly in the time domain to o btain the 
tunnel rates and hence the quasiparticle density,'^^'^ but 
in a functioning detector it is more practical to use the 
time-averaged phase shift of the RF carrier reflected by 
the oscillator to measure the quasiparticle density and 
detect an incoming photon. 

In a parallel-element tank circuit with inductance L 
and capacitance C coupled to a transmission line through 
a capacitance Cc, as shown in Fig. 1, the phase of a 
reflected wave is given by tan^ = —2\Z\Zq/ {\Z\^ — Zq) 
where 
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is the impedance of the tank circuit and Zq ~ 50 ^1 is the 
impedance of the transmission line. In this expression, 
u is the angular drive frequency, uq = 1/ \JL{C + Cq) is 
the center frequency, and cjc = 1/ ^/LCc- Modulation of 
the phase shift with the quantum capacitance comes in 
through modulation of cjq in Eq. (|4|. 

The quantum capacitance of the SCB is proportional 
to the second derivative of the SCB energy with respect 
to the gate charge. When the SCB is tuned to its degen- 
eracy point in the even state, the change in capacitance 
when a single quasiparticle tunnels is given by 



Cl AEc 
Cy\ Ej 



f4Ec 
V Ej 



-3/2N 



C^g ^Ec 
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for an SCB in the charge-eigenstate limit Ec ^ Ej. In 
this expression, Cg is the coupling capacitance between 
the oscillator and the SCB island, while Cy: is the total 
capacitance of the SCB. Ideally, one would like to design 
the SCB and oscillator parameters so that the phase shift 
6(j) induced by SCq is 180°. In the test device discussed 
below, where the oscillator Q ^ 3000, the observed phase 
difference between the even and odd states was approxi- 
mately S<p = 140°, as can be seen in Fig. 2. 

When performing reflectometry with a high-Q oscil- 
lator, the observed quantity is the time-averaged phase 
shift {(/)). The oscillator itself performs a moving average 
with a window length equal to the inverse of its band- 
width, and further averaging can be performed in soft- 
ware or with an external circuit. At all but the lowest 
temperatures and shortest time scales, the phase shift 
switches stochastically between its even- and odd-state 
values in a homogeneous two-rate Poisson random tele- 
graph process with additive Gaussian noise. Since such a 
process is ergodic, we may treat the observed phase shift 
as the ensemble average 
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where (;ieven,odd are the phase shifts in the even and odd 
SCB parity states, and the deviation of the observed 
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phase shift from its "dark" value is ^ = {(j)) — ^even = 
^odd^0- In this way it is quite straightforward to extract 
the quasiparticle density from the phase shift of the re- 
flected wave. 



C. Dark Counts 

In the development of single Cooper-pair devices for 
quantum computation, quasiparticle tunneling has gen- 
erally been regarded as an unwanted effect which is a 
challenge to reproducibly control. As a result, it is rea- 
sonable to expect that in QCD devices there will be sig- 
niflcant quasiparticle tunneling even when no signal pho- 
tons are present. In applications where it is necessary 
to control the number of dark counts, quasiparticle tun- 
neling can be suppressed by creating an energetic barrier 
for quasiparticle tunneling between the absorber and the 
island. This is typically realized in practice by creating a 
gradient in the superconducting gap energy between the 
absorber and the island, either by controlling impurity 
concentrations^^ or fllm thicknesses.^^ However, the 
use of such techniques are likely to limit the sensitivity of 
the detector, requiring a careful engineering tradeoff. Of 
course, it is also possible that the number of dark counts 
can be reduced signiflcantly through improved control 
over materials and more careful electromagnetic flltering 
and shielding. 



III. NOISE SOURCES AND DETECTOR 
SENSITIVITY 

A. Sensitivity 

Given that single-charge devices such as the SCB are 
exquisitely sensitive to the presence of non-equilibrium 
quasiparticles, one of the key advantages of the QCD 
scheme is its extreme sensitivity to electromagnetic ra- 
diation at submillimeter wavelengths. For a material 
with superconducting gap energy A, an incident photon 
with energy hv will generate a number of quasiparticles 
^qp = T\hvl IS. in the absorber, where the factor « 0.57 
is the efliciency with which the energy of the initial photo- 
electron is downconverted into quasiparticles.^^ Assum- 
ing that the distribution of quasiparticles in the absorber 
is spatially uniform (that is, that the dimensions of the 
absorber are small compared to the quasiparticle diffu- 
sion length) this will produce a non-equilibrium quasi- 
particle density nqp = N^^j^ where ^ is the absorber 
volume. 

The sensitivity of the QCD system can be quantifled 
by the noise-equivalent power (NEP), which is deflned as 
the radiant power per square root bandwidth required to 
achieve a signal-to- noise ratio (SNR) of unity at a given 
signal modulation frequency. Since the flnal signal in the 
QCD is the mean phase shift in an oscillator, all relevant 



sources of noise can be expressed as phase noise. For a 
phase noise spectral density S(^{uj\ the NEP is given by 
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where Tqp is the quasiparticle lifetime, Pg is the incident 
radiant source power, and is the ringdown time of the 
oscillator. In this equation, the responsivity 
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is simply deflned as the overall average phase shift per 
quasiparticle, as in the MKID.^ 

The source power is related to the quasiparticle density 
through the rate equation 
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where R is the recombination constant. In aluminum, 
R = 9.6 /im^/s.^The steady-state solution is given by 



(10) 



The phase noise in the QCD signal arises from a variety 
of physical mechanisms. The contributions to the NEP 
from the dominant sources of noise are discussed below. 



B. Telegraph Noise 

In the QCD scheme, the quasiparticle density in the 
absorber is sampled by measuring the rates of quasiparti- 
cle tunneling in the SCB, which occurs in discrete events. 
In practice, the observed quantities are the occupation 
probabilities in the even and odd parity states, which are 
extracted from the average phase shift of the oscillator. 
Since this average must be performed in time over the 
random telegraph signal intrinsic to quasiparticle tun- 
neling, a signiflcant source of phase noise in the QCD is 
"telegraph noise", or quasiparticle shot noise. 

The telegraph noise in the QCD can be characterized 
as a two-state random process with two rates, corre- 
sponding physically to tunneling from the absorber to the 
island (Fin) and from the island to the absorber (Fout)- 
For the purposes of this paper, we treat the tunneling of 
quasiparticles in both directions as homogeneous Poisson 
processes, with no correlation between tunneling events 
and no cross-correlation between quasiparticles tunneling 
into and out of the SCB. At low temperatures and short 
time scales, this assumption has been observed both theo- 
retically and experimentally to be invalid due to a failure 
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of quasiparticles on the island to reach thermal equilib- 
rium. However, at 100 mK, the proposed operation tem- 
perature of the detector, the time scales of these non- 
Poissonian effects are negligibly short and the approx- 
imation of uncorrelated quasiparticle tunneling is well 
justifiedP 

The noise spectrum of a two-rate random telegraph 
signal can be quickly calculated by computing the Fourier 
transform of the autocorrelation function. It is found 
to be 



where F^] = Fi^ 



rinrout/rs 
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(11) 



"out, Tin = Kuqp and Fout is indepen- 
dent of riqp. This is a Lorentzian-type shot noise spectral 
density, which consists of a flat distribution at low fre- 
quencies with a cutoff at the sum transition rate F^] . Note 
that if Fin = Tout, we recover the textbook Lorentzian 
expression for a single-rate random telegraph process. 



quasiparticle states begin to be activated, leading to in- 
creased tunneling rates which swamp out the response of 
the detector to the non-equilibrium quasiparticles gener- 
ated in photon absorption. In closely related devices, 
T* - 250 mK.^o The fluctuation in the QCD phase 
shift signal due to the thermal excitation and recombina- 
tion of equilibrium quasiparticles is known as generation- 
recombination (GR) noise, with an associated frequency- 
independent NEP™ 



NEPgr = 2A 



-A/2/cbT 



(13) 



where the equilibrium quasiparticle number Nl is defined 
in Sec. Ila. This noise arises simply from the fluctuation 
in quasiparticle density from equilibrium generation and 
recombination events, which are assumed to be uncorre- 
lated. As shown in Fig. 6, GR noise does not impose a 
significant limitation on detector sensitivity at the pro- 
posed QCD operating temperature of 100 mK. 



Fano Noise 



E. Other Sources of Phase Noise 



Another key noise mechanism in the QCD is Fano 
noise, which arises due to the uncertainty in the num- 
ber of quasiparticles generated in the absorber by an in- 
dividual photon. When an incident photon is coupled 
to the antenna, a fast photoelectron is generated in the 
absorber, which rapidly thermalizes via electron-phonon 
collisions into an equilibrium distribution of quasiparti- 
cles with mean quasiparticle number Nqp = riqpVt. Since 
this relaxation process is highly correlated, the quasipar- 
ticle generation statistics are sub-Poissonian. The vari- 
ance in the quasiparticle number is not Nqp but FNqp^ 
where the Fano factor F 0.2 quantifies the degree of 
suppression of the fluctuations.^^ In the QCD, the noise- 
equivalent power due to the Fano effect is given by 



V 



(l+c^V|J(l+c.V2). (12) 



While the Fano effect is typically considered as a lim- 
itation on the energy resolution of a detector, the high 
sensitivity of the SCB to quasiparticle density makes this 
effect a small but important contribution to the overall 
NEP. 



D. Generation-Recombination Noise 

At the target QCD operating temperature of 100 mK, 
equilibrium quasiparticle states are strongly frozen out, 
so equilibrium quasiparticle tunneling is exponentially 
suppressed. As a result, the quasiparticle tunneling dis- 
cussed above is due strictly to non-equilibrium quasipar- 
ticles. At a characteristic temperature T*, equilibrium 



In addition to telegraph noise and fluctuations in the 
quasiparticle density, there are a variety of physical ef- 
fects which result in additional noise in the QCD phase 
shift signal. Both lumped-element LC oscillators and 
coplanar waveguide quarter-wave resonators suffer from 
intrinsic phase noise, which has been found to arise 
from coupling to a Fermionic bath of two-level charge 
fiuctuat ors (T LFs) located on substrate and material 
surfaces .'^^'^ This intrinsic oscillator phase noise is typi- 
cally characterized by a l//-type spectrum. This source 
of noise is currently the limiting factor for the sensitiv- 
ity of MKID devices, and is also a key issue limiting the 
performance of phase qubits and nanomechanical oscilla- 
tors. Ongoing progress in reducing this noise has focused 
on the investigation of novel fabrication techniques and 
materials, as well as improved device design. 

Another closely related source of phase noise in the 
QCD is charge noise in the SCB. This is caused by cou- 
pling of the SCB island to the bath of TLFs in the charge 
degree of freedom. This effectively acts as a fluctuating 
gate charge, which is transformed into a fluctuating phase 
shift in the SCB readout. This fluctuation is minimized 
when the SCB is operated at its charge-phase degeneracy 
point, although in a multiplexed array such operation is 
impractical with a single tuning line. 

Finally, noise in the amplifier and readout electronics 
will contribute additional phase noise to the QCD signal. 
State-of-the-art cryogenic HEMT amplifiers have a noise 
temperature < 1 K at 400-800 MHz or 5 K at 4-8 GHz.l^ 
All of the above sources of phase noise can be measured 
in the QCD as additive noise on top of the telegraph 
signal. 

The overall aggregate phase noise, including contribu- 
tions from oscillator phase noise, SCB charge noise, am- 
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FIG. 4: (Color online) Measured aggregate phase noise in a 
test device oscillator with an input power of -125 dBm. This 
device was fabricated using lower- Tc Al/Ti/Au trilayer leads 
to suppress quasiparticle tunneling. The estimated phase 
noise includes contributions from oscillator phase noise, SCB 
charge noise, and noise in the amplifiers and electronics. The 
blue curve is a periodogram of time-domain phase noise mea- 
sured using an oscilloscope. The black line is a fit to a 1/f 
curve with a prefactor of 6.5 degrees at 1 Hz. 



plifier noise, and noise from mixers and other electronics, 
has been measured in a test SCB device and is shown in 
Fig. 4. In this device, which is not coupled to an an- 
tenna, quasiparticle tunneling is suppressed by forming 
the leads of the device from an Al/Ti/Au trilayer with a 
superconducting transition temperature of 450 mK. This 
creates a large energetic barrier for quasiparticle tunnel- 
ing, suppressing the telegraph noise altogether. This al- 
lows us to measure the phase noise directly in the time 
domain, by recording data with an oscilloscope in 10^- 
point frames and computing the Fourier transform. The 
blue line in Fig. 4 shows a 100-average periodogram of 
this data, taken with a 10 kHz amplifier bandwidth. The 



black line is a fit to y S^^^^^{uj) = a/uj^ where the pref- 
actor a = 6.5 degrees at 1 Hz. This fit is later used with 
Eq. ^ to estimate the phase noise NEP in a QCD device 
with practical parameters. 



IV. PERFORMANCE ESTIMATES 

Figures 5-12 illustrate the theoretical performance of 
the QCD detector using an experimentally feasible set 
of device parameters. In choosing detector parameters, 
engineering tradeoffs must be made between detector sen- 
sitivity, saturation power, and practicality of fabrication. 
The selected parameters were chosen as a reasonable 
compromise between these goals, with an eye toward ap- 
plications in far-infrared and sub-millimeter astronomy. 

Figure 5 shows a plot of the mean oscillator response 
^odd^0 as a function of incident power at the detector. 
This is the overall phase shift (j) in the oscillator, time- 



FIG. 5: (Color online) Theoretical average phase response of 
the SCB readout oscillator as a function of detector loading 
power. Solid (blue) curve: Response with small absorber, Q = 
0.1 /im^, optimized for maximum sensitivity. Dashed (red) 
curve: Response with larger absorber, Q = 10 /xm*^, optimized 
for higher saturation power or shorter wavelength radiation. 



averaged over the intrinsic random telegraph noise, as 
given by Eq. The response is shown for two differ- 
ent absorber volumes, 1] = 0.1 /im^ and Q = 10 /im^. 
Both curves were evaluated at a fixed temperature T = 
100 mK, SCB trap depth 6E/kB = 0.1 K, and tunnel 
barrier conductance Gat = 66 jiiU. We have assumed 
tank circuit parameters such that the phase shift between 
parity states is 6(j) = 180°. We have also assumed the 
material parameters A/ks = 2.1 K, Tqp = 110 /is and 
R = 9.6 jum^/s. This figure illustrates that the range of 
detector operation can be tailored simply by choosing an 
appropriate absorber volume, trading off sensitivity for 
saturation power. Although the response is linear over 
at least one decade in power, the simple functional form 
of Eq. (|3| permits detector operation over a wider range. 

Figure 6 shows the theoretical sensitivity of the de- 
tector as a function of signal modulation frequency for 
the design parameters used to compute Fig. 5. The ab- 
sorber volume used was Q = 0.1 /im^, which is the de- 
sign optimized for maximum sensitivity. The incident 
power was taken to be Pg = 10~^^ W. The plot illus- 
trates NEP due to the four dominant noise mechanisms 
described in section III as a function of signal modulation 
frequency. Since the noise arising from different physical 
mechanisms is uncorrelated, the solid (black) line is the 
mean-square sum of the individual NEP values for each 
mechanism. As can be seen from Eq. ([7|), the functional 
form of each NEP curve is dominated by the spectral 
density of noise at low frequencies, and is ultimately lim- 
ited at high frequencies by the quasiparticle lifetime and 
the response time of the readout oscillator. 

The dashed (blue) line is the NEP due to telegraph 
noise, as found by substituting Eq. ( 11 ) into Eq. The 
shape of this curve is dominated by the Lorentzian form 
of the noise spectral density, which is flat at low mod- 
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FIG. 6: (Color online) Noise-equivalent power as a func- 
tion of signal modulation frequency for a variety of phys- 
ical noise mechanisms, as described in the text. Simula- 
tions are performed with fixed temperature T — 100 mK 
and incident power = 10-19 W. Solid (black) line: mean- 
square sum of the NEP due to each noise mechanism. Dashed 
(blue) line: NEP due to telegraph noise. Dotted (green) 
line: NEP due to aggregate excess phase noise. Dash-dotted 
(red) line: NEP due to Fano noise. Solid horizontal (gold) 
line: NEP due to generation-recombination noise. 



ulation frequencies and falls off sharply at the effective 
tunneling rate T^. The particular value of Tj^ = 26 kHz 
shown in Fig. 6 emerges from the particular choice of the 
incident power P^, and will change with optical loading 
in keeping with Eq. ([lO| . For this particular set of device 
parameters, Tout = 15 kHz. It is important to note that 
since the observed signal in the QCD is the phase shift in 
the reflected wave averaged over the telegraph noise, it is 
meaningless to operate the detector at a modulation fre- 
quency faster than the cutoff. As a result, operation must 
be restricted to a modulation frequency 1/r^ <C Tout- 
For the remainder of the calculations in this paper, we 
have taken the operation time to be r^n = 10 ms. 

The dotted (green) line is the NEP due to the ag- 
gregate excess phase noise as a function of frequency, 
adapted from the measurements shown in Fig. 4. The 
plot is of Eq. ^ with ^^/S^{uJ) = a/co where a = 6.5 
degrees at 1 Hz. Note that this noise source is domi- 
nant at low modulation frequencies, where the detector 
is likely to be operated. However, the magnitude of this 
noise is not currently at a fundamental limit, and can 
likely be improved by using better amplifiers and elec- 
tronics, along with lower- noise substrates and materials. 
This will significantly improve the sensitivity of the de- 
tector at low frequencies. The dash-dotted (red) line is 
the NEP due to the Fano noise, as estimated by Eq. (12). 



Note that for the small absorber volume considered here, 
the Fano noise is actually equal to or in excess of the 
telegraph noise. Finally, the solid horizontal (gold) line 
is the NEP due to generation-recombination noise, which 
in the approximation of Eq. ( 13 ) is taken to be frequency- 
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FIG. 7: (Color online) Noise-equivalent power as a function 
of operating temperature for a variety of physical noise mech- 
anisms, as described in the text. Simulations are performed 
with fixed measurement time Tm = 10 ms and incident power 
Ps = 10~^^ W. The color scheme and line styles are identical 
to those of Fig. 6. 



independent. At 100 mK, the design operation temper- 
ature of the device, equilibrium quasiparticle tunneling 
is strongly frozen out and the GR noise is a negligible 
contribution to the total. 

Figure 7 shows the theoretically estimated NEP as a 
function of operating temperature, for the same physi- 
cal noise mechanisms discussed in Fig. 6. In this case, 
the operating time = 10 ms, 1^ = 0.1 /im^, and all 
other parameters are the same as in figs. 5 and 6. The 
color scheme for the plots is identical to that of Fig. 6 
for NEP due to telegraph noise, excess phase noise, Fano 
noise, and generation-recombination noise. At very low 
temperatures, it is interesting to observe that the NEP 
increases sharply. This is due to the sharp increase in the 
island-to-absorber tunneling rate Tout which occurs at 
low temperatures, where quasiparticles in the SCB island 
fail to reach thermal equilibrium and tunnel elastically.'^ 
This sharply decreases the responsivity of the detector at 
low temperature, increasing the NEP. Also note that in 
the proposed region of operation, 100 mK, the sensitivity 
of the detector is quite flat. Above 150 mK, tunneling 
of quasiparticles in thermal equilibrium begins to dom- 
inate, and the generation-recombination noise becomes 
the dominant noise source. Note that the Fano noise is 
independent of temperature, since the temperature de- 
pendence enters into the NEP through the responsivity, 
which cancels in the derivation of Eq. (12). 

Figure 8 shows the detector NEP as a function of 
loading power, calculated using the same expressions as 
Figs. 6 and 7. The design parameters are the same as 
those used in Fig. 7, with the temperature held at 100 
mK and Q = 0.1 /im^. The solid (blue) line is the total 
NEP, which is the mean-square sum of the NEPs from 
telegraph noise, excess phase noise, Fano noise, and GR 
noise. Note that the NEP increases dramatically as the 
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FIG. 8: (Color online) Noise-equivalent power as a function 
of loading power, as described in the text. Simulations are 
performed with fixed modulation time Tm = 10 ms and op- 
erating temperature T = 100 mK. Solid (blue) line: mean- 
square total NEP due to all noise mechanisms. Dashed (green) 
line: Photon shot noise at a wavelength A = 30 /im. Dotted 
(red) line: Photon shot noise at A = 1 mm. 



incident power approaches 10~^^ W. At this power level, 
the detector is approaching saturation, as can be seen 
from Fig. 5. 

For far-infrared and submillimeter radiation, another 
key qualification is whether the detector is sensitive 
enough to be limited by the shot noise of the incident 
signal itself. To compare the QCD sensitivity with the 
shot noise limit, we have plotted the NEP due to pho- 
ton shot noise NEP^ = ^J2Pshc|\ at two wavelengths, 
A = 30 /im and A = 1 mm, as shown by the dashed 
(green) and dotted (red) lines, respectively. From this 
figure, one can see that the QCD detector is predicted 
to be shot-noise limited over a wide range of operating 
powers. 

In the same vein. Fig. 9 shows a plot of the saturation 
power of the detector as a function of frequency com- 
pared to the shot noise limits. In this plot, the solid 
(blue) curve is the saturation power -Psat? defined as the 
power at which the response ^ of the detector falls within 
one noise standard deviation a = y/^^{uJ)uo of the phase 
response limit 5(j)^ evaluated at a frequency uo. The over- 
all structure follows the frequency dependence shown in 
Fig. 6, with saturation power decreasing at low frequen- 
cies due to 1 // phase noise at at high frequencies due to 
the finite quasiparticle lifetime and finite response time 
of the readout oscillator. The dashed (green) and dot- 
ted (red) lines are an illustration of the critical power 
above which the detector ceases to be shot-noise limited, 
as illustrated in Fig. 8. The points where the shot noise 
limits fall rapidly indicate frequency ranges over which 
the detector is not shot-noise limited at all. The dashed 
and dotted curves in Fig. 9 correspond to the shot noise 
limits for A = 30 /im and A = 1 mm, respectively. 

As we have stressed above, a critical QCD design pa- 




Modulation Frequency (Hz) 

FIG. 9: (Color online) Saturation power as a function of fre- 
quency, for the sample parameters used in Fig. 8. Solid (blue) 
curve: Physical saturation power of the detector, as defined in 
the text. Dashed (green) curve: Critical power at which the 
detector ceases to be shot noise limited, as defined in the text, 
for A = 30 iim. Dotted (red) curve: Shot noise limit for A = 1 
mm. The points where the shot noise limits drop rapidly off 
scale indicate frequency ranges where the detector is not shot 
noise limited at all. 
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FIG. 10: (Color online) Noise-equivalent power as a func- 
tion of absorber volume Q, for fixed measurement time Tm = 
10 ms, operating temperature T = 100 mK, and incident 
power Ps — 10~^^ W. The color scheme and line styles are 
identical to those of Figs. 6 and 7. While the sensitivity of the 
detector improves with a smaller absorber, engineering trade- 
offs must be made with saturation power and practicality of 
fabrication. 

rameter is the absorber volume Vl. Fig. 10 shows theo- 
retical NEP as a function of ^ for the four noise sources 
considered above. The colors and line styles for the plot 
are the same as for Figs. 6 and 7. As can be seen from 
Fig. 5, the sensitivity of the device improves as the ab- 
sorber volume is reduced, assuming that the absorber- 
island tunneling rate depends linearly on the quasiparti- 
cle density. This improvement in sensitivity clearly comes 
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FIG. 11: (Color online) Comparison of QCD detector sensitiv- 
ity as a function of wavelength with the detector requirements 
for the BLISS spectrometer. Detector requirements and spec- 
trometer loading data are adapted from Ref. 15. 



at the expense of saturation power, so Q must be tai- 
lored for a specific application. The absorber cannot be 
made arbitrarily small, since it must be coupled to the 
antenna and must be comparable in volume to the SCB 
island. Furthermore, the absorber must be large enough 
for the generated quasiparticle population to equilibrate 
below the gap edge of the Nb antenna within the dif- 
fusion time, to ensure quasiparticle confinement in the 
absorber region. In the calculations discussed above, we 
have set O = 10~^^ m^, which is a reasonable compromise 
between these concerns. As with the temperature depen- 
dence, note that the Fano noise is volume-independent. 



V. APPLICATIONS 

Such sensitive submillimeter radiation detectors are 
critically needed for future experiments in far-infrared as- 
trophysics. Approximately half of the total light emitted 
from stars and black hole accretion over the history of the 
universe has been absorbed by dust and reradiated in the 
10-1000 jam band. However, detailed spectroscopic inves- 
tigations of this spectral region have remained difficult 
due to the challenges of ground-based observation and a 
need for sensitive detectors. While significant advances 
have recently been made with germanium photoconduc- 
tors and beam-isolated silicon nitride TES bolometers, 
the QCD could play an important role as a sensitive de- 
tector which can be frequency-multiplexed into large ar- 
rays in a straightforward fashion. 

A proposed instrument which is a natural fit for 
the QCD detector is the Background-Limited Infrared- 
Submillimeter Spectrometer (BLISS), a spectrograph de- 
signed for the Space Infrared Telescope for Cosmology 
and Astrophysics (SPICA).^^ SPICA is a 3.5 m, 4.5 K 
cold space telescope planned for launch in 2013. One of 



the key technological challenges for realizing the BLISS 
spectrometer is the development of a sensitive submil- 
limeter detector array. In Fig. 12, we show the detec- 
tor sensitivity requirements for the BLISS spectrome- 
ter alongside the theoretical NEP for the QCD detector, 
plotted as a function of incident radiation wavelength. 
For each wavelength, we have taken the design optical 
loading power for the BLISS spectrometer and computed 
the theoretical NEP as in Fig. 8, using the QCD design 
parameters given above. For most of the BLISS spectral 
region, the QCD detector meets these sensitivity require- 
ments by a wide margin. The sharp increase in NEP at 
long wavelengths is due to increased optical loading in 
the spectrometer design. At A = 30 /im, the design load- 
ing power is Pg = 1.62 x 10~^^ W, while at A = 800 /im, 
the loading power Pg = 5.4 x 10~^^ W. 

Because the QCD sensitivity can be tuned to the quasi- 
particle density by varying the absorber volume, the 
QCD may also be useful in applications at shorter wave- 
lengths. Furthermore, there exist a wide variety of non- 
astrophysical applications for sensitive detectors in far- 
infrared and submillimeter regime, including earth sci- 
ence, planetary science, biomedical technology and de- 
fense applications. 



VI. CONCLUSIONS 

We describe a scheme for a sensitive pair-breaking 
photon detector applicable for far-infrared and sub- 
millimeter radiation. Incident photons coupled with an 
antenna break Cooper-pairs in a superconducting ab- 
sorber. A single Cooper-pair box in contact with the 
absorber is used as a probe to measure the density of 
quasiparticles. The SCB signal can be read out by using 
a radio-frequency capacitance measurement or by em- 
bedding the SCB in a microwave resonator. This scheme 
lends itself naturally to large-scale frequency multiplex- 
ing. The sensitivity and performance of the proposed 
device has been calculated, including several physical 
noise sources. These sources include telegraph noise, 
fano noise, generation-recombination noise, and excess 
resonator, charge and electronic phase noise. The detec- 
tor has very favorable predicted sensitivities, with a min- 
imum NEP on the order of 10~^^ W/V^Hz, and can be 
readily fabricated using existing techniques. The QCD 
promises excellent sensitivity, and is a natural tool for 
applications in far-infrared astrophysics. 
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